Comment on "Central limit behavior in deterministic dynamical systems" 



OO 

o 
o 

(N 

Oh: 

CD 
00 



o 

CD 

e 

I 



O 
o 



> 

O 

o 

00 

o 



X 



Peter Grassberger 

Complexity Science Group, Department of Physics and Astronomy, University of Calgary, Calgary, Canada 

(Dated: September 8, 2008) 

We check claims for a generalized central limit theorem holding at the Feigenbaum (infinite 
bifurcation) point of the logistic map, made recently by U. Tirnakli, C. Beck, and C. Tsallis (Phys. 
Rev. 75, 040106(R) (2007)). We show that there is no obvious way that these claims can be made 
consistent with high statistics simulations. We also refute more recent claims by the same authors 
that extend the claims made in the above reference. 

PACS numbers: 05.20.-y,05.45.Ac,05.45.Pq 



In [H, the authors claim that it is "less known in the 
physics community" that the central hmit (CL) theorem 
holds for deterministic but mixing dynamical systems. 
Our first comment is that this connection is perfectly 
well known and was never doubted. An early reference 
is i. 

But my main criticism concerns the treatment of the 
(non-mixing) logistic map at the Feigenbaum (infinite bi- 
furcation 2]) point in [l|. Let us consider trajectories 
Xi+i ^ fa{xi) of length A^, with fa{x) = a - . The 
Feigenbaum point is at a = Ac = 1.40115518909 . . .. Fol- 
lowing [H], we study sums Y — 'Y^^=n^+i''"^ ^^'^ their 
distributions for random at a = and for large N . 
Here, Nq is the length of a possible discarded transient. 
At first we shall consider the case Aq = 0, i.e. no tran- 
sient is discarded. We understand that this is also the 
case studied in [l[. At least the authors do not men- 
tion any transient, although there are reasons (discussed 
later) to suspect that that they might have used A^o > 0. 

Denoting by (y) the average over ccq, the claim of [3| 
is that the centered and suitably rescaled sums 
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for 7 = 1.5, with b w 4/3 and c ~ 0.1 [g]. Moreover, it 
is claimed that the same distribution, with identical 7, 6, 
and c, is found for the modified logistic map fa.z(,x^ — 
a — with z = 1.75 and z = 3. If true, this universality 
would be remarkable. 

Unfortunately, none of the above claims seem to be 
correct. It is straightforward to do the necessary simu- 
lations to estimate p{y). In all simulations, xo was uni- 
formly distributed in [0, a]. Results, for several values of 
a at and slightly above Cc, are shown in Fig.l. Indeed, 
in this figure are shown histograms of the non-rescaled 
and non-shifted sums Y, for A^ — 16384. Similar re- 
sults were obtained for other values of A^. They have 
markedly different behavior left and right of the central 
peak Yc « 86333. For F < 1^ we observe a very steep 
rise, P{Y) ~ e°'*^^, while the decrease for Y > Yc is 
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1: Histograms of Y for N = 16384 and for four values 
at or slightly above criticality. 



much more gentle, P{Y) '-^ e 



Superposed on both 



exponentials are periodicities which obviously result from 
the hierarchical structure of the Feigenbaum attractor. 

In obtaining Fig.l we have not discarded any transient 
(A'o = 0). Notice that the entire dynamics is transient 
at the Feigenbaum point [1, |3| . Thus discarding a finite 
transient would introduce a new time scale and ruin any 
hope for scaling. 

This would be different, if values a > Oc were consid- 
ered as in and, presumably, also in [l|: The results 
of Ref . ^ are most easily understood as artifacts gener- 
ated by using the 9-digit approximation 1.40115519 for 
ttc- Let us consider a value of a where the attractor con- 
sists of n = 2*^ "bands" 0, Orbits on it jump peri- 
odically between the bands, but are chaotic within each 
band. Then it makes sense to discard transients of length 

(n) 

^ n. On each band, the n— fold iterated map fa is mix- 
ing. Thus y is a sum over n series of random variables, 
each of which shows normal CL behavior for N ^ 00. 
Therefore, Y also shows normal CL behavior in the limit 
N — > 00, n = const. 

Deviations from normal CL behavior can be expected 
only when taking a joint limit a — > Oc (i.e., n — > 00) and 
A^ ^ 00. Indeed, in Q it was proposed that N ^ v? . 
Notice, however, that in this case we are no longer dealing 
with the problem posed in the central limit theorem, i.e. 
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FIG. 2: Distributions of yl\J var{y), normalized to p{0) = 1, 
for various values of A'^ and n, where a is set to the n ^ n — 1 
band merging point. In all cases, A^o > 16384. Similar results 
were obtained also for other values of a. Notice that the 
statistics in any curve of this figure is at least 10 times higher 
than in any of the curves in [ij, |5| 
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FIG. 3: Distributions of y/y^ var{y) for = 65536, normal- 
ized to p(0) = 1, plotted against y^ + const on a log- log plot. 
As in Fig. 2, long transients (A^o ~ 65536 in most cases) have 
been discarded, and the control parameter a of the logistic 
map is chosen as the n ^ n/2 band merging point. Accord- 
ing to Eq. (2), one would expect straight lines with slopes 
—b. Apart from the rather unsystematic deviations at large y 
which could be effects which vanish in the limit N ^ n —> oo, 
one sees a systematic downward curvature for intermediate 
y and strong systematic upward (downward) curvatures for 
n > \/7V (n < y/N) at very small y. Notice that one has two 
curves for each n, one for y > and one for y < 0. 

the asymptotics of partial sums in an infinite sequence 
of random variables. Thus, strictly speaking, we are no 
longer deahng with (normal or abnormal) central limit 



behavior at all. 

In the following we shall, for definiteness, only deal 
with band-merging points, where n = 2'^' bands merge 
into n/2 bands as a is increased. But similar behavior is 
found also for other values of a. Indeed, when looking at 
distributions of y for large n, N » n, and iVo >> n, 
one finds heavy-tailed distributions, see Fig. 2. But as 
closer inspection shows, they are in general not described 
by Eq. (2) (see Fig. 3) . Apart from the steps and discon- 
tinuities at large y which might recede to infinity in the 
limit indicated above, the main deviations are: 

• A systematic downward curvature in Fig. 3 for in- 
termediate to large y; 

• Deviations from straight line behavior at very small 
y, both for rt << \/N and for n >> a//V. If at aU, 
the data are compatible for small y with Eq. (2) 
only for a very narrow region of Njn}. 

In view of this it seems very unlikely that the rough agree- 
ment with Eq. (2) is more than a numerical coincidence. 
The data shown in Figs. 4 and 5 of 1] are definitly not 
well fitted by Eq. (2) (not even for very small y/cr). We 
might add that equally good (or bad) fits would be ob- 
tained with Levy stable distributions [7| , which moreover 
have more theoretical justification. 
Some final remarks: 

• The behavior described here is seen only when iV is 
a power of 2. Otherwise, one observes completely 
different behavior. 

• The fluctuations of Y arc, for Nq,N ^ n ^ 1, 
tiny. All structures shown in Figs. 2 and 3 (includ- 
ing the tails!) extend, before centering and dividing 
by Nf, over a range AY < 10'^. For N = 65536, 
this is to be compared to (Y) « 34533, i.e. all rela- 
tive fluctuations are smaller than 3 x 10^* [8]. The 
reason for this is that the motion on an n— band 
attractor with large n is extremely regular, with 
the chaos conflned to very narrow bands. Thus if 
a generalized CL theorem holds for this problem 
in any sense, it is completely unobservable in any 
experimental situation. 

• Since the phenomenon illustrated in Figs. 2, 3 
seems to describe corrections to the scaling limit 
of the Feigenbaum map, it is not clear how much it 
depends on the original map one starts from and on 
the distribution of xq. The only phenomenon dis- 
cussed in this comment which has a realistic chance 
to be experimentally accessible and is likely to be 
universal is the behavior shown in Fig. 1. It is dom- 
inated by chaotic transients, and is very far from 
anything described in Refs. [ll, 
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